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ABSTRACT 


Production  planning  for  large-scale  production 
systems  requiring  the  allocation  of  numerous  re¬ 
sources  is  considered.  It  is  demonstrated  how  the 
dynamic  activity  analysis  developed  by  Shephard 
leads  to  linear  programming  solutions  of  produc¬ 
tion  planning  problems.  Three  types  of  planning 
problems  are  formulated:  maximization  of  output 
levels  for  a  given  time  horizon;  minimization  of 
production  duration  for  given  output  histories; 
and  minimization  of  production  costs  for  given  out¬ 
put  histories. 


Production  Planning  for  Multi-Resource  Network  Systems 


Robert  C.  Leachman 
University  of  California,  Berkeley 


1.  Introduction 

Previous  efforts  (von  Lanzenauer  [4]  and  Candea  [1])  to  mathematically  model  capaci¬ 
tated,  multistage  production  systems  have  been  motivated  by  manufacturing  shop  environ¬ 
ments,  in  which  many  products  are  to  be  produced  using  a  given  network  of  facilities.*  The 
problem  considered  is  to  determine  workforce  levels  and  product  lot  sizes  in  each  time  period 
so  as  to  minimize  costs  to  meet  external  demand  schedules  [I]  [4], 

The  focus  of  this  paper  is  planning  for  production  systems  in  which  the  production  net¬ 
work  elements  are  dedicated  to  producing  a  single  product,  but  allocation  of  numerous 
resources  among  the  elements  is  required,  and  other  kinds  of  production  planning  problems  are 
posed.  Shephard  et.  al.  [7]  (81  have  developed  a  continuous  flow  dynamic  activity  analysis 
model  of  production,  in  which  a  network  of  activities  characterizes  the  component  tasks  of  pro¬ 
duction.  Required  facilities  and  other  resources  are  considered  as  inputs  to  be  allocated  among 
the  activities.  The  presentation  in  [8]  is  taken  as  an  appropriate  point  of  departure  here.  In  the 
following,  this  model  is  extended  to  include  inventory  capacities  as  in  [4],  initial  inventories  of 
intermediate  products,  and  classes  of  exogenous  inputs. 

Three  types  of  production  planning  problems  are  formulated  and  solved  using  linear  pro¬ 
gramming  methods.  The  problem  types  considered  are  maximization  of  output  accumulations 
by  a  given  horizon;  minimization  of  production  duration  for  required  output  histories;  and 
minimization  of  costs  for  given  output  requirements. 


*  No  attempt  is  made  here  to  reference  multistafe  mode  line  efforts  of  (incapacitated  systems  or  pure  serial  or 
parallel  structures,  for  a  survey  of  such  « (forts,  see  Candea  III. 
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2.  The  Motel 

Following  Shephtrd  and  Al-Ayat  [81.  the  production  system  is  viewed  as  a  network  of 

production  activities  which  are  denoted  by  Ax,  A2 . A*.  In  the  network,  nodes  represent 

activities  and  arcs  indicate  intermediate  product  transfers,  ie,  the  use  of  each  activity’s  output  as 
input  by  other  activities.  The  operation  of  each  activity  A,  is  measured  in  terms  of  an  intensity 
Junction  z,U).  r— 0, 1.2.  •  •  • .  whose  value  at  time  t  indicates  activity  input  during  fr.r-4-1)  and 
output  at  time  (r+1)  when  taken  with  technical  coefficients  defined  as  follows: 

(a)  c,(/),  r-0, 1,2 . i— 1 . N.  where  c,(t)  is  the  amount  of  output  of  activity  A,  at 

time  /  per  unit  intensity  of  activity  A,. 

(b)  a,k(t),  /— 0. 1,2 . *-1.2 . NK,  /-I . N,  where  «*(/)  is  the  amount  of  exo¬ 

genous  input  type  *  required  at  time  t  per  unit  intensity  of  activity  A,.  The  first  NS  <  NK 
inputs  are  designated  non-storable  resources  which  cannot  be  accumulated;  the  remaining 
exogenous  inputs  can  be  accumulated,  and  are  termed  storable  resources. 

(c)  5,y(f),  /— 0.1.2 . /—I . N.  y— 1 . N,  where  a„U)  is  the  amount  of  inter¬ 

mediate  product  from  activity  A,  required  at  time  t  per  unit  intensity  of  operating  A,. 

For  a  time  horizon  T  for  production  activity,  we  introduce  the  following  technical  limita¬ 
tions  on  the  system: 

I  w -0.1.2 . r-i 

?,(/>J  .the  activity  intensity  bounds,  natural  bounds  resulting  from  available 

workspace  and  other  limitations  not  considered  as  exogenous  input; 

It »— o.  i .  2 . r-i 

.  the  time  histories  of  non-storable  resource  levels  available  for 

Ja-i . vs 

input  to  the  system; 

(»  r— 0. 1 . r-l 

Yk  (t)  J  .the  time  histories  of  storable  resources  made  available  to  the  sys- 

Ja-vs+i . v* 

tern,  where 

t  Mr> 

r-O 

is  the  cumulative  amount  of  resource  *  supplied  during  (O./l; 
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(d)  |//»v,°|  ,  the  initial  inventories  of  activity  product  for  intermediate  uses;  and 

11  i-o.  i . r-i 

cap,  0)1  ,  the  bounds  on  accumulations  of  activity  product  awaiting  inter* 

J—i . .v 

mediate  uses,  arising  from  limited  storage  capacity  for  in-process  inventories. 

For  our  purposes  here,  the  intensity  z,(f)  of  activity  A,.  /-I . N,  on  each  time  inter¬ 
val  1/./+1),  /- 0, 1,2 . r-1.  shall  be  partitioned  into  effort  producing  intermediate  product, 

z/(f),  and  effort  producing  final  product,  t ) .  where 

z/(f)  +  z,f(f)  -  z,(f).  (1) 

These  variables  indicate  the  allocation  of  activity  output  produced  during  [f,f— 1)  to  final  and 
intermediate  uses. 

A  production  plan  is  a  specification  over  some  finite  period  [0.  T)  of  the  activity  intensities 

1i-o . r-i 

. 

i-l . v 

Such  a  plan  is  said  to  be  feasible  for  UT)  if  the  plan  belongs  to  the  set  L(T)  defined  by  the  fol¬ 
lowing  inequalities: 

L(T)1.  £ a,t<f)[z/(f)  +  z'(/)l  <  **(/),  *-l . NS.  f-0,1 . T-I. 

—i  1  ' 

L(T)2.  t,  £<i,*(t)[z/(t)  +  z,'(t)|  <  £  Kt(r>.  *-/VS+ 1 . NK. 

r-0  .-I  1  1  r-0 

. r-i. 

L(T)3.  ±  ia/l(T)(z,,(r)  +  z,f(r)l  -  £c,(t+I)z/(t)  <  invf. 

»-l  ,-l  1  1  r-0 

>-l . N.  f-l . f-1.  and 

£a„(0)\zH0)  A-  z,f(0)l  <  invf,  j-\ . N. 

/-i  1  1 

L(T)4.  'jjcj(r+l)z/(r)  -  tth  (tllz^T)  +  zf(r)|  <  cap  fit)  —  invf. 


t—O  r-l  i-l 

y-1 . N.  f-l . r-l.  and 
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-  £a,,(0)|z/(0)  +  z,f(0)  1  <  capj( 0)  -  invf,  j—  1 . M 

<-i  1  1 

L(T)5.  z/(t)  +  zfU)  <  Ij(t).  j- 1 . jV,  /- 0 . r-1. 

zjU ),  z/(/)  >  0.  y-1 . M  /-0 . r-i. 

Constraints  L(T)1  and  L(T)2  express  resource  limitations.  Constraints  L(T)3  insure  ade¬ 
quate  intermediate  product  transfers  occur  to  support  production  activity,  while  constraints 
L(T)4  insure  that  inventories  of  intermediate  products  do  not  exceed  capacities.  Finally,  con¬ 
straints  L(T)S  limit  intensities  to  non-negative  values  less  than  intensity  bounds. 

The  set  of  linear  inequalities  L(T)  constitutes  a  continuous  flow  model  of  production,  in 
which  any  positive  intensity  of  activity  operation  supplies  completed  product  to  final  or  inter¬ 
mediate  uses,  or  to  inventory.  In  the  case  that  intermediate  products  of  a  system  are  large, 
discrete  units,  precedence  relations  occur  between  activities  output  unit  by  output  unit,  and 
constraints  L(T)3  and  L(T)4  must  be  modified.  This  case  will  not  be  treated  here,  and  the 
reader  is  referred  to  16],  in  which  a  dynamic  activity  analysis  was  developed  on  a  critical  path 
analysis  network. 

In  the  case  more  than  one  activity  produces  a  certain  product,  constraints  LIT) 3  and 
L(T)4  must  be  modified  for  the  activities  in  question.  See  [6]  and  [8].  However,  with  such 
revisions  LIT)  still  constitutes  a  set  of  linear  inequalities.  For  simplicity  of  exposition,  we 
assume  in  what  follows  that  no  two  activities  produce  the  same  product. 

3.  Production  Programming 

3.1.  Output  Maximisation 

In  this  section,  programs  are  formulated  for  the  maximization  of  value  or  mix  functions 
of  final  output  accumulations.  We  consider  first  the  case  where  a  specific  product  mix  of  final 
output  is  desired,  and  the  problem  is  to  maximize  the  scale  of  this  mix  accumulated  by  a  time 
horizon  T. 

Let  zv+,  be  a  variable  indexing  the  scale  of  the  accumulation  by  time  T.  The  amounts  of 
the  various  products  will  be  related  by  coefficients 

fli„v+t.  <“l. . . . .  N, 


where  a,  v+1zv+i  is  the  amount  of  final  product  from  activity  A,  accumulated.  The  problem  in 
question  is  formulated  as  a  linear  program  as  follows. 

Maximize  zv+i 

subject  to 

Ol.  5,..v+i*v+i  <  2  c,U+l)z,FU).  /— 1 . N. 

(-0 

Ir-o . r-i 

€  L(T). 

<-i . ,v 

03.  Z\+i  ^  0. 

In  general,  the  program  involves  [S(N)  +  NKHT)  +  N  +  1  variables,  and 
[3(/V)  +  NK](T)  +  M  constraints.  Clearly,  the  time  horizon  (ie,  the  number  of  time  periods) 
is  the  most  sensitive  factor  in  terms  of  problem  size  which  can  be  handled.  The  structure  of  the 
constraint  set  can  be  modified  by  rewriting  constraints  L(T)3  and  L(T)4  in  terms  of  intermedi¬ 
ate  product  inventory  variables 

|«iv/(r)J.  r-0 . r-I, 

which  are  the  slack  variables  defined  by  constraints  L(T)3.  Using  these  variables,  we  reformu¬ 
late  constraints  L(T)3  and  L(T)4  as  follows: 

L(T)3.  Ja;((r)[z/(f)  +  z,f(f)|  -  c/(r)z/(r — 1)  —  invfit- 1)  +  invj(t)  -  0, 

/-i  1  * 

7  -1 . N,  r— 1 . r-1.  and 

£fl/((0)|z/(0)  +  zf  (0)  1  +  mv/(0)  -  mv®,  7-1 . N. 

/-1  1  1 

L(T)4.  invj(t)  <  cap,lt),  y-1 . N.  r-0 . r-1. 

With  this  revision,  it  is  evident  that  the  constraints  L(T)1,  L(T)4  and  L(T)S  apply  only 
time  period  by  time  period,  and  the  constraint  matrix  exhibits  partial  block  diagonal  structure. 
Potential  is  thus  offered  for  application  of  large-scale  programming  procedures  such  as 
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TT 


decomposition.  (See  (5].) 

We  next  consider  the  case  where  the  vaiue  of  output  produced  is  to  be  maximized.  We 
suppose  each  product  /  has  a  constant  unit  price  p,.  The  maximum  value  of  output  accumulated 
from  production  activity  during  (O.D  is  then  given  by  the  optimum  of  the  following  linear  pro¬ 
gram. 

r-i  s 

Maximize  £  £  P,c,(t+l)z^(t) 

0  /“I 

subject  to 

(l'-o . r-i 

zf(t),  z'(r)|  tUT). 

The  remarks  about  problem  size  and  structure  concerning  the  previous  program  apply  here  as 
well,  as  only  the  product  mix  variable  and  the  N  constraints  Ol  have  been  deleted. 

3.2.  Time  Minimization 

In  this  kind  of  planning  problem,  there  are  final  output  demands  which  must  be  met,  but 
the  overall  production  duration  is  to  be  minimized.  Final  output  demands  are  expressed  in 
cumulative  terms  as  follows.  Let 

£/,(/),  f— 1,2 . r,  1-1 . /V, 

denote  the  required  cumulative  delivery  of  final  product  /  by  time  t.  Here,  we  are  considering 
the  situation  where  early  delivery  of  final  products  is  acceptable  or  even  desirable.  These 
demands  imply  constraints 

t  c,(T+l)z(f(T)  >  U,U),  /- 1 . N.  r— 1 . r.  (2) 

r*0 

We  first  consider  the  problem  of  finding  the  latest  starting  time  for  production  activity 
sufficient  to  satisfy  (2).  A  feasible  production  plan  for  this  problem  would  satisfy  the  linear  ine¬ 
qualities  (2)  and  L(T).  An  optimal  plan  would  have  the  characteristic  that 

z,(0)  -  r,(l)  -  z,(2)  -  •••  -  z,(,0)  -  0.  i-1 . N.  (3) 


where  r0  >s  as  large  as  possible.  Such  a  plan  may  be  found  (if  one  exists)  by  solving  a  sequence 
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of  Phase  I  linear  programs  (see  i2])  as  follows. 

The  set  of  inequalities  under  consideration  is  of  the  form 

A*  +  Bx  -  b 
t  >  0.  x  >  0. 

where 


i-  |z{  (0),  zf<0) . z((0),  z«0>.  z(  (1 ),  zf(l) . zC(l).  z((l). 

z{  ( 7* — 1).  zf(7-l) . z((  7— 1).  Zv(T— 1)  J; 


. xm+n)  ; 


m-(N)(T).  n-(NK)(T)  +  3(<V)(7).  /-2(JV)<7)  ; 


(4) 


A  is  the  (ffi+n)x/ matrix  of  coefficients  of  activity  intensities  in  (2)  and  LIT),  where  the  first 
m  rows  arise  from  (2); 

B  is  the  (m+n)x(m+n)  matrix  of  coefficients  of  slack  variables  for  said  constraints;  and 
b  is  the  (m+n)  vector  of  right  hand  side  constant  terms  of  the  constraints. 

For  the  inequalities  organized  in  this  fashion,  the  solution  algorithm  is  presented  below: 
Step  0.  Initialize  t  -  7-1. 

Step  1.  Solve  the  Phase  I  problem  with  the  first  (/ V)(t )  columns  of  the  tableau  corresponding 
to  (4)  deleted.  If  a  feasible  solution  is  found,  stop;  then 

I  \'-o . r- 1 


z/O).  z,f(r) 


f-i. 


is  an  optimal  production  plan.  If  the  problem  is  infeasible,  go  to  Step  2. 

Step  2.  If  r-0,  stop;  then  the  set  of  inequalities  is  infeasible.  Otherwise,  decrease  r  to  t-1  and 
go  to  Step  1 . 

The  algorithm  is  seen  to  initially  ignore  all  columns  in  the  tableau  corresponding  to 
activity  intensities  in  periods  before  time  (7-1),  and  to  then  attempt  to  find  a  basic  solution.  If 
none  can  be  found,  columns  corresponding  to 

z/(  7-2),  z,f(7- 2).  /-I . N. 


are  also  considered.  The  algorithm  continues  to  allow  the  use  of  columns  corresponding  to 
activity  operation  one  time  period  before  the  earliest  period  of  activity  operations  allowed  by 
the  previous  iteration.  The  algorithm  terminates  either  the  first  time  a  feasible  basic  solution  is 
found,  or  else  all  columns  are  adjoined  without  finding  one.  In  the  former  case,  an  optimal  pro¬ 
duction  plan  is  found,  and  in  the  latter  case,  the  output  schedule  (2)  is  infeasible  for  the  limita¬ 
tions  L(T). 

We  next  consider  the  problem  of  finding  the  earliest  time  all  product  accumulations  can 
be  completed.  It  is  immediately  apparent  that  an  approach  similar  to  that  considered  above  can 
be  used  to  solve  this  problem.  A  sequence  of  Phase  I  procedures  is  again  suggested,  but  in  this 
case  starting  with  the  possibility  of  positive  activity  intensities  only  during  [0. 1),  and  proceed¬ 
ing  forwards  in  time.  Later  production  activity  is  allowed  period  by  period  until  either  the  first 
time  a  feasible  basic  solution  is  found,  or  else  the  horizon  is  reached  without  finding  one.  In 
the  former  case,  an  optimal  production  plan  is  obtained,  and  in  the  latter  case,  the  output 
schedule  is  infeasible  for  the  limitations  imposed. 

3.3.  Cost  Minimization 

In  this  section  we  formulate  the  problem  of  determining  a  minimum  cost  production  plan 
which  meets  a  given  final  output  schedule  expressed  in  the  form  of  (2). 

Non-storable  resources  are  assumed  to  have  capacity  costs  corresponding  to  the  peak 
demands  for  each  such  resource.  These  resources  cannot  be  accumulated,  so  that  the  produc¬ 
tion  system  must  have  the  capability  to  accommodate  peak  loads.  Storable  resources,  however, 
have  prices;  these  resources  account  for  the  variable  cost  of  production  activity.  It  is  assumed 
that  storable  resources  can  be  procured  as  required,  so  that  inventories  of  same  are  ignored.  We 
assume  intermediate  product  inventories  also  have  capacity  costs,  corresponding  to  peak  storage 
requirements.  These  inventories  will  also  bear  holding  costs  representing  opportunity  charges 
for  unproductive  capital. 

Assuming  linear  capacity  costs,  the  problem  is  formulated  as  a  linear  program  as  follows. 
Let 

Cv(/)  -  (cf  (r) . CVr) )• 


be  the  vector  of  costs  per  unit  capacity  for  non-storable  resources  maintained  during  (r,r+l);  let 
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CO) -(cW'> . C*</>). 

be  the  price  vector  for  storable  resources  procured  for  use  during  lr.r+1);  let 

C -(f)  -  (c-J  (r) . Cv(/)|. 

be  the  vector  of  costs  per  unit  storage  capacity  maintained  during  (f.t+1)  for  intermediate  pro¬ 
ducts;  and  let 

HO)  -  (//,</) . //vO)). 

be  the  holding  costs  for  intermediate  products  held  during  [/,r+D. 

To  serve  as  variables  in  the  minimization,  let 

x  -  (*, . *4 

denote  the  peak  requirements  in  any  unit  time  interval  of  non-storable  resources;  let 

cap  -  I  cap, . capX 


denote  the  required  intermediate  product  inventory  storage  capacities;  and  let 

1 1-0 . T- 1 


wiv/O).  r/( t ) ,  r,0 1 ) 

be  the  inventory  and  intensity  variables  as  before. 
For  given  intensity  bounds 

I  f»0. 


2,0) 


i-l. 


. .  r-i 

.  v 


and  initial  intermediate  product  inventories 


the  minimum  cost  production  plan  meeting  the  output  schedule  (2)  is  given  by  the  optimum  of 
the  linear  program 


Minimize  £C*(/)-X+£  f  £  Q(.t)alkU)\z/U)  +  zf(t)\ 

/- o  /-o  *-vs+ 1  <-i  1  1 

+  X‘  C-(/)cap  +  £  H,U)inv!U) 

r«-0  /— 0  » ■  I 

subject  to 

Cl.  £  c,(T+l)z,f(r)  >  £/,</).  /-I . Af,  r-1 . T. 

r-0 

C2.  £  fl,*(/)[r/0)  +  r,f(/)l  -  Xk  <  0,  *-l . NS.  r-O, 1 . f-1. 

,-i  1  J 

C3.  X0„(/)[r/(r)  +  2'^^1  ~  cl(l)zj(l—\)  -  invj(l—\)  +  invjU)  —  0, 

<-i  1  1 

7- 1 . AT.  /-I . r-1.  am/ 

Xa„(0)[r/(0)  +  zfCO)]  +  mv/(0)  -  invf.  >-l . M 

,-i  1  1 

C4.  «•>/(/)  -  cap,(r)  <  0.  /-I . A',  r-0 . r-1. 

C5.  z((t)  +  zf(t)  <  ?,(/).  7-1 . N.  r- 0 . r-1. 

C6.  X  -  [*, . JTvsj  >  0. 

cap  -  jcflP! . copy)  ^  0. 

inv'O)  —  | inv{(t) . /tayv(r)j  ^  0,  r— 0 . T—  1. 

r/(f),  ;/"(/)  ^  0.  7-1 . M  t-0 . r-1. 

Here  constraints  C2  define  the  required  non-storable  resource  capacities,  and  constraints 
C4  define  the  required  intermediate  product  storage  capacities.  Constraints  C3  and  CS  deal  with 
inventory  balance  and  intensity  bounds  in  the  same  manner  as  the  treatment  of  output  maximi¬ 
zation  problems,  while  constraints  Cl  repeat  (2). 

In  general,  the  program  includes  (47V  +  NS)(T+\)  variables  and  (4N  +  NS)(T)  con¬ 
straints.  As  before,  the  fineness  of  the  time  grid  is  the  most  sensitive  factor  in  terms  of  the 
problem  size  which  can  in  practice  be  solved.  A  bordered  angular  configuration  for  the 
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constraint  matrix  is  now  displayed,  in  which  constraints  C2,  C4,  and  CS  exhibit  a  block  diagonal 
structure  with  coupling  variables.  Although  this  is  a  more  difficult  structure  than  that  for  the 
output  maximization  problems,  nonetheless  it  can  be  exploited.  See  [Si. 

As  an  alternative  to  the  constant  capacities  for  each  non-storable  resource  defined  by  con¬ 
straints  C2,  one  may  allow  capacities  to  be  adjusted  from  time  period  to  time  period  according 
to  linear  costs.  Many  authors  have  formulated  labor  workforce  levels  in  this  fashion,  allowing 
hiring  and  firing  in  each  period.  See  for  example  [3]  or  [4],  Such  formulations  may  be 
integrated  here  as  appropriate. 
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